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ABSTRACT 

We present the exact analytic expressions to compute, assuming the emitted Gamma-Ray Burst (GRB) ra- 
diation is not spherically symmetric but is confined into a narrow jet, the value of the detector arrival time at 
which we start to "see" the sides of the jet, both in the fully radiative and adiabatic regimes. We obtain this 
result using our exact analytic expressions for the EQuiTemporal Surfaces (EQTSs) in GRB afterglows. We 
re-examine the validity of three different approximate formulas currently adopted for the adiabatic regime in 
the GRB literature. We also present an empirical fit of the numerical solutions of the exact equations, compared 
and contrasted with the three above approximate formulas. The extent of the differences is such as to require a 
reassessment on the existence and entity of beaming in the cases considered in the current literature, as well as 
on its consequences on the GRB energetics. 

Subject headings: gamma rays: bursts — ISM: kinematics and dynamics — relativity 



1. INTRODUCTION 



After the work by 'M ao & Yil jl994 . the possibility that 
Gamma-Ray Bursts (GRBs) originate from a beamed emis- 
sion has been one of the most debated issue about the nature of 
the GRB sources in the current literature (see e.g. Piran^2QII4; 
iMeszaros 2006, and references therein). I n particular , on th e 
ground of the theoretical considerations bv' SarietalJ l fT999l) . 
it was conjectured that, within the framework of a conical jet 
model, one may find that the gamma-ray energy releas ed in all 
GRBs is narrowly clustered around 5 x 10^" ergs iFrail et alJ 
12001 . ^ ^ 

In a recent letter JBianco & Ruffinil200 5bll) we analyzed the 
approximate power-law relations between the Lorentz gamma 
factor and the radial coordinate usually adopted in the current 
GRB literature. We pointed out how such relations are found 
to be mathematically correct but only approximately valid in 
a very limited range of the physical and astrophysical param- 
eters and in an asymptotic regime which is reached only for 
a very short time, if any. Therefore, such relations were there 
shown to be not applicable to GRBs. Instead, the exact ana- 
lytic solutions of the equations of motion of a relativistic thin 
and uniform shell expanding in the interstellar medium (ISM) 
in the fully radiative and adiabatic regimes were there pre- 
sented. 

This program of identifying the exact analytic solutions in- 
stead of approximate power-law solution is in this letter car- 
ried one step forward. Using the above exact solutions, we 
here introduce the exact analytic expressions of the relations 
between the detector arrival time tf^ of the GRB afterglow 
radiation and the corresponding half-opening angle d of the 
expandi ng source visible are a due to the relativistic beaming 
(see e.g. lRuffini et al.l2003 al). Such visible area must be com- 
puted not over the spherical surface of the shell, but over the 
EQuiTemporal Surface (EQTS) of detector arrival time f^, i.e. 
over the surface locus of points which are source of the ra- 
diation reac hing the observer at the same arrival time (see 
iBianco & Ruffini. .2004 .2005a. for details). The exact ana- 
lytic expressions for the EQTSs in GRB afterglows, which 
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have been presented in 'Bianco & Ruffini' ("2005a"), are there- 
fore crucial in our present derivation. This approach clearly 
differs from the ones in the current literature, which usually 
neglect the contributions of the radiation emitted from the en- 
tire EQTS. 

The analytic relations between and d presented in this 
Letter allow to compute, assuming that the expanding shell 
is not spherically symmetric but is confined into a narrow jet 
with half-opening angle do, the value {t'^)jei of the detector 
arrival time at which we start to "see" the sides of the jet. A 
corresponding "break" in the observed light curve should oc- 
cur later than (?„);« (see e.g. Sari et al. 1999). In the current 
literature, (?„ ))e, is usually defined as the detector arrival time 
at which 7 ~ 1 /i9o, where 7 is the Lorentz factor of the ex- 
panding shell (see e.g. Sari et al. 1999, and also our Eq.(|3 
below). In our formulation we do not consider effects of lat- 
eral spreadings of the jet. 

In the current literature, in the case of adiabatic regime, 
different approximate power-law relations between {t^)jet 
and do hav e been pr esented, in contrast to each other (see 
e.g. Sari et al. 1999; Panai tescu & M eszaros 1999: Pan aitescul 
2006). We show here that in four specific cases of GRBs, en- 
compassing more than 5 orders of magnitude in energy and 
more than 2 orders of magnitude in ISM density, both the one 
by Panaitescu & Meszaros ( 1999) and the one by Sari et ak] 
( 1999) overestimate the exact analytic result. A third relation 
just presented by Panaitescu ( 2006) slightly underestimate the 
exact analytic result. We also present an empirical fit of the 
numerical solutions of the exact equations for the adiabatic 
regime, compared and contrasted with the three above approx- 
imate relations. In the fully radiative regime, and therefore in 
the general case, no simple power-law relation of the kind 
found in the adiabatic regime can be established and the gen- 
eral approach we have outlined has to be followed. 

Although evidence for spherically symmet ric emission in 
GRBs is emerging from observations (Soderberg et al. 200' 
and from theoretical argumentations (Ruffini etal. 2004; 
120051) . it is appropriate to develop here an exact theoretical 
treatment of the relation between {tg)jet and do- This will 
allow to make an assessment on the existence and, in the pos- 
itive case, on the extent of beaming in GRBs, which in turn is 
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going to be essential for establishing their correct energetics. 

2. ANALYTIC FORMULAS FOR THE BEAMING ANGLE 

The boundary of the visible region of a relativistic thin and 
uniform shell expanding in the ISM is defined by (see e.g. 
iRuffini et al..2003a. and references therein); 



cos-i? = - , 
c 



(1) 



where -d is the angle between the line of sight and the radial 
expansion velocity of a point on the shell surface, v is the ve- 
locity of the expanding shell and c is the speed of light. To find 
the value of the half-opening beaming angle -do correspond- 
ing to an observed arrival time it^)je,, this equation must be 
solved together wijji the equation d escribin g the EQTS of ar- 
rival time (tfi)jet (jBianco & Ruffinil l2005aft . In other words, 
we must solve the following system: 



cosi9o = ^ 

COSTCO = COS 



\ ■ (2) 

'J<^'1 \EQTS[(tf,)j„]f 

It should be noted that, in the limit i?o — + and v c, this 
definition of (t^)jei is equivalent to the one usually adopted in 
the current literature (see sec. [Q. 

2.1. The fully radiative regime 

In this case, the analytic solution of the equations of motion 
gives (see Bianco & Ruffini 2005a b): 
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1 + (M,-.,„,/Mb) + (Mis,,, /Mb 



1 + {Mism/MB) (1 +7o') [1 + 5 (M,,,„/Mb)] 



(3) 



where 70 and are respectively the values of the Lorentz 
gamma factor and of the mass of the accelerated baryons at 
the beginning of the afterglow phase and M,j„, is the value of 
the ISM matter swept up to radius r. Mism = (47r/3)OTpn,j„,(r-'- 
/•q-'), where ro is the starting radius of the baryonic matter 
shell, nip is the proton mass and n,j„, is the ISM number den- 
sity. Using the analytic expression for the EQTS given in 
iBianco & RuffiniU2005a^ ■ Eq.© takes the form: 
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where to is the value of the time f at the beginning of the 
afterglow phase, m° = {4/3)TTmpnis^rl, r* is the initial size of 
the expanding source, A = [(M^ -/«°)/m°]'/^, C = Ms^i'jo- 
1 )/ (70 + 1 ) and z is the cosmological redshift of the source. 

2.2. The adiabatic regime 

In this case, the analytic solu tion of the equations of motion 
gives (see .Bianco & Ruffini 2005albh : 



7o-l 7o + 



Mb 



(5) 




10"^ 10^ 10^ 
Detector arrival time {ta)jg, (s) 

Fig. 1. — Comparison between the numerical solution of Eq.|4| assum- 
ing fully radiative regime (blue line) and the con'esponding one of Eq.O 
assuming adiabatic regime (red line). The depaiture from power-law behav- 
ior at small arrival time follows from the constant Lorentz 7 factor regime, 
while the one at large angles follows from the approach to t he non-relativistic 
regime (see details in section|4]and Fig. |4] as well as in IBianco' &Ruffiml 
l2005hh . 



Using the analytic exp ression for the EQTS given in 
IBianco & Ruffini! J2005al) . Eq.(|3 takes the form: 



COSt?o 



COSl^o 




(6) 
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Kl+z) r ^7^-1 

where all the quantities have the same definition as in Eq.©. 

2.3. The comparison between the two solutions 

In Fig. [2 we plot the numerical solutions of both Eq.Q, 
corresponding to the fully radiative regime, and Eq.(|6j, corre- 
sponding to the adiabatic one. Both curves have been plotted 
assuming th e same initial conditi ons, namely the ones of GRB 
991216 rsee lRuffini et alJ2003ah . 

3. COMPARISON WITH THE EXISTING LITERATURE 

Three different approximate formulas for the relation be- 
tween (f^)je/ and do have been given in the current literature, 
all assuming the adiabatic regime. iPanaitescu & MeszarosI 
in"999.) proposed: 



cosi?o=i 1-5.9 X 10' 
ISari et al instead, advanced 

do lA X 
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. 1+2 _ 



3/8 



(7) 



(8) 



In both Eq.(0 and Eq.(|8j, (tf,)jet is measured in seconds, E is 
the source initial energy measured in ergs and n,jm is the ISM 
number density in particles/cm^. The formula bv iSari et alJ 
( 119991) has been applied quite ofte n in the curren t literature 
(see e .g. IFrail et ah 120011: .Ghii'landa etalJ 12004 iFox et alJ 
I200I. 
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Fig. 2. — Comparison between the numerical solution of Eq.|6| (red line) and the corresponding approximate formulas given in Eq.JsJ (blue line), in Eq.0 
(black line), and in Eq.JgJ (green line). All four curves have been plotted for four different GRBs: a) G RB 991216 (see Ruffini et al. 2003a), b) GRB 980519 
(see lRuffini et al. 2006). c) GRB 031203 (see Bemardini et al. 2005), d) GRB 980425 (see Rufiini et all20d4h) . The ranges of the two axes have been chosen to 
focus on the sole domains of application of the approximate treatments in the current literature. 



Both Eq.Q and Eq.(|8j compute the arrival time of the pho- 
tons at the detector assuming that all the radiation is emitted 
at ■& = (i.e. on the line of sight), neglecting the full shape 
of the EQTSs. Recen tly, a new expression has been proposed 
by Panaitescu (200&), again neglecting the full shape of the 
EQTSs but assuming that all the radiation is emitted from 
I? = 1/7, i.e. from the boundary of the visible region. Such an 
expression is: 



:5.4x 10' 









. 1 + z . 



3/8 



(9) 



In Fig. |2]we plot Eq.Q, Eq.(|8} and Eq.(|9} together with the 
numerical solution of Eq.(|6j relative to the adiabatic regime. 
All four curves have been plotted assuming the same initial 
conditions for four different GRBs, encompassing more than 
5 orders of magnitude in energy and more than 2 orders of 
magnitude in ISM density: a) GRB 991216 ( see Ruffi ni et aP 
l2Q03a) . b) GRB 980519 (see Ru ffini et al.h2006) . c) GRB 
031203 (se e Bernardini et al. 2003), d) GRB 98 0425 (see 
iRuffini et al . 2004b). The approximate Eq.(|8j bv ISari et all 
il999t) and Eg. by Panaitescu (2006) both imply a power- 
law relation between -do and {t^)jei with constant index 3/8 
for any value of -do, while Eq.(0 bv .Panaitescu & MeszarosI 
lll999h implies a power-law relation with constant index 3/8 



only for -do ^ (for greater -do values the relation is trigono- 
metric). 

All the above three approximate treatments are based on the 
approximate power-law solutions of the GRB afterjglow dy- 
namics which have been shown in lBianco & Ruffinil (I2005bl) 
to be not applicable to GRBs. They also do not take fully 
into account the structure of the EQTSs, although in different 
ways. Both Eq.Q and Eq.(|S), which assume all the radiation 
coming from = 0, overestimate the behavior of the exact 
solution. On the other hand, Eq.(|9}, which assumes all the ra- 
diation coming from i? ~ I/7, is a better approximation than 
the previous two, but still slightly underestimates the exact 
solution. 

4. AN EMPIRICAL FIT OF THE NUMERICAL SOLUTION 

For completeness, we now fit our exact solution with a suit- 
able explicit functional form in the four cases considered in 
Fig. 12] We chose the same functional form of Eq.(|9}, which is 
the closer one to the numerical solution, using the numerical 
factor in front of it (i.e. 5.4 x lO-') as the fitting parameter We 
find that the following approximate expression: 



^0=^5.84x10^^ 



3/8 



(10) 
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Analytic solution in adiabatic regime ■ 
Fitting tunction - 




Detector arnval time (tgjjg, (s) 

Fig. 3. — The overlapping between the numerical solution ot 'Eq . i^ (thick 
green lines) and the approximate fitting function given in Eg. 1101 (thin red 
lines) is shown in the four cases (a-d) represented in Fig. |2] 




lo" 10^ 10^ 
Detector arrival time {ta)jg| (s) 

Fig. 4. — Comparison between the numerical solution of Eq.iSl (think 
green lines) and the approximate fitting function given in Eg. 1101 (thin red 
lines) in all the four cases (a-d) represented in Fig. |2] The ranges of the two 
axes have been chosen to have their full extension (i.e. the one of Fig. fTl . 
The dashed gray lines are the boundaries of the region where the empirical 
fitting function can be applied. 



is in agreement with the numerical solution in all the four 
cases presented in Fig. |2](see Fig. |3}- However, if we en- 
large the axis ranges to their full extension (i.e. the one of 
Fig. 0, we see that such approximate empirical fitting for- 
mula can only be applied for -do < 25° and {ti)jet > 10^ s (see 
the gray dashed rectangle in Fig. 0}. 



An equivalent empirical fit in the fully radiative regime 
is not possible. In this case, indeed, there is a domain in 
the ({tg)jet,i}o) plane where the numerical solution shows a 
power-law dependence on time, with an index ^ 0.423 (see 
Fig. [0. However, the dependence on the energy cannot be 
factorized out with a simple power-law. Therefore, in the fully 
radiativ e regime, which is th e relevant one for our GRB model 
(see e.g. lRuffini et al.l2003a t). the application of the full Eq.@ 
does not appear to be avoidable. 



5. CONCLUSIONS 



We have presented in Eqs.( l4l6> the exact analytic relations 
between the jet half-opening angle i?o and the detector ar- 
rival time (t^)jei at which we start to "see" the sides of the 
jet, which may be used in GRB sources in which an achro- 
matic light curve break is observed. The limiting cases of 
fully radiative and adiabatic regimes have been outlined. Such 
relations differs from the approximate ones presented in the 
current hterature in t he adiabatic regime: both the ones by 
iPan aitescu & Meszaros ( .1,999) and bv .Sari_et__aL. (199£) over- 
estimate the exact analytic result, while the one just presented 
b V iPanaitescuii (i2006 ) slightly underestimate it. 

For a limited domain in the ((f^') je, ,-do) plane defined in Fig. 
13 and only in the adiabatic regime, an empirical fit of the nu- 
merical solution of the exact Eq.(|6j has been given in Ea.(llO>. 
However, in the fully radiative regime such a simple empirical 
power-law fit does not exist and the application of the exact 
Eq.(|3 is needed. This same situation is expected also to occur 
in the general case. 

In light of the above results, the assertion that the gamma- 
ray energy re leased in all GRB s is narrowly clustered around 
5 X 10^" ergs (i Frail et al.l20 01') should be reconsidered. In ad- 
dition, the high quaUty data by Swift, going without gaps from 
the "prompt emission" all the way to latest afterglow phases, 
will help in uniquely identifying the equations of motion of 
the GRB sources and the emission regimes. Consequently, 
on the ground of the results presented in this Letter, which 
encompass the different dynamical and emission regimes in 
GRB afterglow, an assessment on the existence and, in the 
positive case, on the extent of beaming in GRBs will be pos- 
sible. This is a step in the determination of their energetics. 



We thank an anonymous referee for his/her interesting sug- 
gestions. 
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